Gradient estimate of a Dirichlet eigenfunction on 
a compact manifold with boundary 



Yiqian Shi* and Bin Xu* f 



Abstract. Let ex (x) be an eigenfunction with respect to the Dirichlet Laplacian An on a compact 
Riemannian manifold N with boundary: An ex = A 2 ex in the interior of N and ex — on the 
boundary of N . We show the following gradient estimate of ex : for every A > 1 , there holds 
A||ex||oo/C < ||Vex||oo < CA||ex||oo) where C is a positive constant depending only on N. In the 
proof, we use a basic geometrical property of nodal sets of eigenfunctions and elliptic apriori 
estimates. 
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1 Introduction 

Let (N, g) be an n-dimensional compact smooth Riemannian manifold with smooth 
boundary 3N and An the positive Dirichlet Laplacian on N. Let L 2 (N) be the space 
of square integrable functions on N with respect to the Riemannian density dv(N) = 



x/g(x) dx := Y / det~TgijT dx. Let ei(x), ez[x) } ■■■ be a complete orthonormal basis in L 2 (N) 
for the Dirichlet eigenfunctions of An such that < A 2 < < A 2 < • • • for the corresponding 
eigenvalues, where ej(x) (j = 1,2, ...) are real valued smooth function on N and Aj are 
positive numbers. Also, let ej denote the projection of L 2 (N) onto the 1-dimensional 
space Cej. Thus , an L 2 function f can be written as f = ^|^o e j (fl> where the partial sum 
converges in the L 2 norm. Let A be a positive real number > 1 . We define the spectral 
function and the unit band spectral projection operator xa as follows: 



e(x,y,A) := Y_ e jM e j(y)> 

Aj <A 

XA f:= Y. e i (f) • 

Aj£(A,A+1] 
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Grieser [5] and Sogge [13] proved the L°° estimate of xa, 

IbCAflloo < CA^^Hflb 



(1) 



where ||f|| T (1 < r < oo) means the L r norm of the function f on N. In the whole of this paper 
C denotes a positive constant which depends only on N and may take different values at 
different places, if there is no otherwise stated. The idea of Grieser and Sogge is to use the 
standard wave kernel method outside a boundary layer of width CA _1 and a maximum 
principle argument inside that layer. By using the maximum principle argument and the 
estimate (1), Xu [15] proved the gradient estimate of xa 

l|V X Aflloo<CA( n+1 '/ 2 ||f|| 2 . (2) 

Here V is the Levi-Civita connection on N. In particular, Vf = Y.) 3f/3xj is the gradient 
vector field of a C 1 function f, the square of whose length equals Y.ij 9 l '(9f/9x0(9f/9xj)- 
One of his motivation is to prove the Hormander multiplier theorem on compact manifolds 
with boundary. Seeger and Sogge [11] firstly proved that theorem by using the parametrix 
of the wave kernel on manifolds without boundary. All the results mentioned in the in- 
troduction have their analog on compact manifolds without boundary. See the details in 
the introduction of Shi-Xu [12] and the references therein. In general, the method used 
in manifolds without boundary is not valid for the problems on manifolds with boundary. 
In particular, on manifolds with boundary the Hormander multiplier theorem cannot be 
obtained by the standard pseudo-differential operator calculus as done on manifolds with- 
out boundary, since the square root of the Dirichlet Laplaican is not a pseudo-differential 
operator any more and one cannot obtain the L°° bounds for xa and Vxa only by using 
the Hadamard parametrix of the wave kernel. 

In the paper, by rescaling xxf at the scale of A -1 both outside and inside the boundary 
layer of width CA _1 , we obtain by elliptic apriori estimates a slightly stronger estimate 
than (2) as follows: 

Theorem 1.1. Let f be a square integrable function on N. Then, for every A > 1, 
there holds 

II V X Af Hoc < C (A||xAf Hoc + A- 1 ||A N XAf Hoc) • (3) 

Remark 1.1. Putting f(-) = Y-^ety a+i] e jM e j(0 i n (1)> we obtain the uniform estimate 
of eigenfunctions for all x e N , 

Y_ lejMI^CA- 1 . 

AjG(A,A+1] 
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Actually it can be proved that those two estimates are equivalent. By using the Cauchy- 
Scahwarz inequality, we obtain the gradient estimate (2) from Theorem 1 together with 
the above inequality. Similarly, that estimate is equivalent to the uniform estimate for all 

xeN, 

Y_ IVejMf^CA*- 1 . (4) 

Aj e(A, A+1] 

Remark 1.2. By the finite propagation speed of the wave equation, the asymptotic 
formula of derivatives of the spectral function e(x,y,A) in Theorem 1 [14] which is proved 
by the standard wave kernel method, also holds for each interior point x of N. Much more 
general asymptotic formulae are given in Theorems 1.8.5 and 1.8.7 of Safarov-Vassiliev 
[10]. In particular, we have the following asymptotic formula that as A — > oo 

^ n+2 - '-T1+1- 



Tl A ' 

L' v qM' 2 = 2(4^r(2T¥) +0 ' (A 



) ) 



where the constant in the reminder term O x (A n+1 ) depends on the distance of x to the 
boundary of N. Hence, the exponents of A in estimates (2) and (4) are sharp at x as 
A — > oo. For each point z on the boundary of N, Ozawa [9] used the heat kernel method 
to show the asymptotic formula that as A — > oo, 

A n+2 



L 

Aj<A 



+ o(A^ +2 ), 



where y is the unit outer normal vector field on the boundary of N. Hence, (2) and (4) 
are also sharp on the boundary. 



Corollary 1.1. Let e\{x) be an eigen] 'unction with respect to the positive Dirichlet 
Laplacian A^ on N: A^e^ = A 2 e A in the interior of N and e\ = on the boundary of 
N . Then, for every A > 1 , there holds the upper bound estimate of Ve A : 

llVeAHoo < CA||e A || 

oo- 



Proof. Putting f = e\ in the estimate (3), we obtain the corollary. □ 

Actually, by the basic geometric property of nodal sets of an eigenfunction, we can 
find a complete picture for the L°° norm of Ve\ in the following: 

Theorem 1.2. Let e\[x) be an eigenfunction with respect to the Dirichlet Laplacian 
operator An on N without boundary: An ex = A 2 ex- Then, for every A > 1 , there holds 

A||e A ||oo/C < ||Ve A ||oo < CA||e A ||oo. (5) 
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Remark 1.3. The authors [12] proved the analog of Theorems 1 and 2 on compact man- 
ifolds without boundary. The proof in this paper is more complicated than there because 
we need to do analysis at points near the boundary. We believe that there also hold the 
analog for k-covariant derivatives V k XA^ and V k e;\ on N. We plan to discuss this question 
in a future paper. 

Remark 1.4. Let i[>j be the normal derivative of e, at the boundary 3N of N. The lower 
bound estimate 

l|4>j||L°°(aN) > C||ej|| L oo( N) 

does not hold in general. Using Examples 3-5 in Hassell-Tao [6] and doing a little bit 
more computations, we can see that the above estimate does not hold on the flat cylinder, 
the hemisphere and the spherical cylinder. We hope to find a sufficient condition for the 
lower bound estimate in a future work. 

We conclude the introduction by explaining the organization of this paper. In Section 
2, we show the lower bound of the gradient VeA by the basic geometrical property of the 
nodal set of eigenfunctions. In Section 3 we use the rescaling method and the Holder 
estimate about elliptic PDEs to show (3) and the upper bound part of (5). The point is 
to do the rescaling both outside and inside the boundary layer of width CA _1 . 

2 Lower bound of VeA 

The nodal set of an eigenfunction e\ of An is the zero set 

Z eA : ={xGN:e A (x)=0}. 

A connected component of the open set N\Z eA is called a nodal domain of the eigenfunc- 
tion eA- We have the same definition for manifolds without boundary. 

Lemma 2.1. (Briining [1]) Let M be a compact Riemannian manifold without bound- 
ary. Let A > 1 and e A be an eigenfunction of the positive Laplacican Am-' Am ^a = A 2 eA- 
Then there exists a constant C only depending on M such that each geodesic ball of 
radius C/A in M must intersect the nodal set Z eA of e A . 

Proof. A proof written in English is given by Zelditch in pp. 579-580 of [16]. □ 
We need a manifold-with-boundary version of Lemma 1 as follows: 
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Lemma 2.2. Let A > and A 2 be greater than the smallest eigenvalue A 2 of the 
Dirichlet Laplacian An. Let e-\ be an eigenfunction of An.' An ex =A 2 eA in the inte- 
rior Int(N) of N and e-\ = on the boundary 3N o/N. Then there exists a positive 
constant D only depending on N such that each geodesic ball of radius D/A contained 
in Int(N) must intersect the nodal set Z eA of e\. 

Proof. We here adapt the proof of Zelditch [16] with a slight modification. 

Step 1 We show the following fact: There exists a constant C such that for each 
interior point p of N and each positive number r > satisfying that the distance d(p, 3N) < 
r from p to 3N is less than r, the smallest Dirichlet eigenvalue A 2 (B(p, r)) of the geodesic 
ball B(x, r) is bounded from above by C/r 2 . Since d(p, 3N) > r, we may assume that 
there exists a geodesic normal coordinate chart (xi,...,x n ) on the ball B(p,r). Let g be 
the Riemannian metric of N on B(p,r) with coefficients = g(3/3xi, 3/3xj) and go be 

the Euclidean metric dx 2 H hdx 2 on B(x, r). Take < Ci < 1 depending only on N so 

that the Euclidean ball B(p,cir; go) is contained in the metric ball B(p,r; g) of N. Since 
Ci < 1 , by the definition of Rayleigh quotient, 

A 2 (B(p,r;g)) <A 2 (B(p, Cl r;g)). 

Since N is compact, by comparing Rayleigh quotients, there exists Ci > depending only 
on N such that 

A 2 (B(p,cir; g)) < c 2 A 2 (B(p,cir; g )). 

On the other hand, by change of variables, we have 

,2,„, ^ A 2 (B(p,1;g )) 
Af (B(p, Cl r; g )) = ^—^ '-. 

Combining the above three inequalities and setting C = C2C 1 ~ 2 A 2 (B(p, 1 ; go)) , we complete 
the proof. 

Step 2 Take a geodesic ball B(p,r) in N such that d(p, 3N) > r. Suppose that it is 
disjoint from the nodal set Z eA . Then it is completely contained in a nodal domain Dj 
of e A . But A 2 = A 2 (Dj) < A 2 (B(p,r)) < C/r 2 . Hence, r < \/CA. Taking D = 2y/C, we 
complete the proof. □ 

Proof of the lower bound part of Theorem 2 Take a point x in N such that 
I^aMI = ll e A||oo- By the Dirichlet boundary condition, the distance d from x to 3N is 
positive. 



5 



Case 1 Assume d > D/A. Then there exists point y in the geodesic ball B(x, D/A) 
with center x and radius D/A such that e A (y) =0. We may assume A so large that there 
exists a geodesic normal chart (r,0) G [0, D/A] x § n_1 (1) in the ball B(x, D/A). By the 
mean value theorem, there exists a point z on the geodesic segment connecting x and y 
such that 



3e A 
3r 



A , . 



D 



l e A| 



Case 2 Assume d < D/A. We may assume A so large that there exists a unique 
geodesic y : [0, d] — > N of arc length parameter connecting x and 3N, 

Y(0)=x, y(d)e3N. 



Since eA(y(d)) = 0, by the mean value theorem, there exists to in (0, d) such that 



de A (y(t)) 



dt 



(to) 



k\MI a 

d -D l|eA||o °- 



□ 



3 Estimate for VxAf 

3.1 Outside the boundary layer 

Recall the principle: On a small scale comparable to the wavelength 1 /A, the eigen- 
function e A behaves like a harmonic function. It was developed in H. Donnelly and 
C. Fefferman [2] [3] and N. S. Nadirashvili [8] and was used extensively there. Recently 
Mangoubi [7] applied this principle to studying the geometry of nodal domains of eigen- 
functions. In this section, for a square integrable function f on N we give a modification 
of this principle, which can be applied to the Poisson equation 

A N XAf= Y- A ? e ' (f) in Int(N) 

AjG(A,A+l] 

with the Dirichlet boundary condition x A f = on 3N. In particular, in this subsection, we 
do the analysis outside the boundary layer L^ A = {z G N : d(z, 3N) < 1/A} of width 1 /A. 

Take an arbitrary point p with d(p, 3N) > 1/A. We may assume that 1/A is suffi- 
ciently small such that there exists a geodesic normal coordinate chart (xi,...,x n ) on the 
geodesic ball B(p, 2/A) in N. In this chart, we may identify the ball B(p, 2/A) with the n- 
dimensional Euclidean ball E(2/A) centered at the origin 0, and think of the function XAf 
in B(p, 2/A) as a function in B(2/A). Our aim in this subsection is to show the inequality 

l(Vxxf)(p)l < C (A||XAf|l Loo(B(2/A)) +A- 1 ||A NX Af|l Loo(B(2/A]) ) • (6) 
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For simplicity of notions, we rewrite u = x^f and v = A N xxf in what follows. The Poisson 
equation satisfied by u in B(1/A) can be written as 

— 7=H (g ij \/g3 Xj u) =v. 
Vg y 

Consider the rescaled functions u\{y) = u(y/A) and V\[y) = v(y/A) in the ball B(2). 
The above estimate we want to prove is equivalent to its rescaled version 

|(Vu A )(0)| < c (||uaII lco(b(2)) +a- 2 KII loo(b(2)) ) • (7) 

On the other hand, the rescaled version of the Poisson equation has the expression, 

Y. 9 vt (9a VjPhjUa) = - a_2 \/9a v a, (8) 

y 

where g iit \[y) = gtjC-y/A), g%[y) = g ij (y/A) and ^/g7(y) = (y^Ky/A)- 

For each < a < 1 , there exists K > such that the C a norm of the coefficients $y/g~\, 
yfgx in B(2) are bounded uniformly from above by K, and the smallest eigenvalue of the 
n x n matrix (g^-y/glOij in B(2) bounds from below by 1/K, for all A > 1. By Theorem 
8.32 in page 210 of Gilbarg-Trudinger [4], there exists constant C = C(n, a, K) such that 

II Ua IIc'.«(b(u) - c (iI Ua IIl-(b(2)) +a_2 H Va IIl-(b(2)))' 

This is stronger than the estimate (7). Therefore, we complete the proof of Theorem 1.1 
outside the boundary layer L^. 

3.2 Inside the boundary layer 

Using the notions in subsection 3.1, We are going to prove the following estimate: 

l|Vu|| L oo( Ll/A ) < C (a||u||oo + A -1 ||v||oo) , (9) 

with which combining (6) completes the proof of Theorem 1.1. 

We may assume that A is sufficiently large so that there exists a geodesic normal 
coordinate chart (z ; , Zn) on the boundary layer L 3 / A ={pgN: d(p, 3N) < 3/A} with respect 
to the boundary 3N. Hence, for each point (z ; , Zn) G L 3 / A , we have < z n < 3/A and 

d((z , ,z n ),3N) =Zn. 

For each point q G 3N and r > 0, denote by B + (q, r) the set of points of N with distance 
less than r to q. Denote by B + (r) the upper half Euclidean ball 

{x = (xi,...,x n ) G R n : |x| < r, x n > 0} 
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centered at the origin and with radius r. Then, for each q £ 3N, there exists a geodesic nor- 
mal chart on B + (q, 3/A) such that the exponential map exp q at q gives a diffeomorphism 
from B+(3/A) onto B+(q,3/A). 

Since {B + (q, 2/A) : q £ 3N} forms an open cover of L^, the question can be reduced to 
showing the analog of (9) on B + (q, 2/A) for each q. We only need to prove its equivalent 
rescaled version, 

ii Vua IIl-(b +( 2)) ^ c (iI Ua IIl-(b + (3)) +a_2 II Va IIl-(b +( 3)))' ( 10 ) 

where ua and va are the the rescaling function of u an v, respectively. Observe that ua 
and va satisfy the Poisson equation (8) in the upper half Euclidean ball B(3) and the 
Dirichlet boundary condition, 

ua=0 on the portion {x £ B(3) : x n = 0} 

of the boundary of B + (3). For each < a < 1 , there exists K > such that the C a norm 
of the coefficients g^'^/gA, y/9x in B + (3) are bounded uniformly from above by K, and the 
smallest eigenvalue of the nxn matrix [g^i/gx^ in B + (3) bounds from below by 1/K, 
for all A > 1. By Theorem 8.36 in page 212 of Gilbarg-Trudinger [4], there exists constant 
C = C(n, a, K) such that 

II Ua IIc'.«(b +( 2)) ^ c (iI Ua IIl-(b + (3)) +a_2 II Va IIl-(b +( 3)))- 

This is a stronger estimate than (10). 
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